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ABSTRACT 

On large angular scales (> 60°), the two-point angular correlation function of the temperature 
of the cosmic microwave background (CMB), as measured (outside of the plane of the Galaxy) by 
the Wilkinson Microwave Anisotropy Probe, shows significantly lower large-angle correlations than 
expected from the standard inflationary cosmological model. Furthermore, when derived from the 
full CMB sky, the two lowest cosmologically interesting multipoles, the quadrupolc [i = 2) and 
the octopole (£ = 3), are unexpectedly aligned with each other. Using randomly generated full-sky 
and cut-sky maps, we investigate whether these anomalies are correlated at a statistically significant 
level. We conclusively demonstrate that, assuming Gaussian random and statistically isotropic CMB 
anisotropics, there is no statistically significant correlation between the missing power on large angular 
scales in the CMB and the alignment of the I = 2 and I = 3 multipoles. The chance to measure the 
sky with both such a lack of large-angle correlation and such an alignment of the low multipoles is 
thus quantified to be below 10 -6 . 



1. INTRODUCTION 

Several anomalies in the cosmic microw ave background 
maps observed by the W MAP satellite (Bennett et al. 
|2003[[Jarosik et aT 2010 ) have been recently observed and 
much discussed in the literature. These include align- 
ments of the lowest modes of CMB anisotropy with each 
other, and with geometry and direction of motion of the 
Solar System, as well as unusually low power at these 
largest scales. Attempts to explain these anomalies in- 
clude astrophysical, instrumental, cosmological causes, 
as well as arguments that faulty data analysis or a pos- 
teriori statistics are at work (for a r eview, see Copi, 
Huterer, Schwarz, & Starkman (2010)). Two particu- 



lariy puzzling features at large ang ular scales are the 
alignments of large-angle anisotropies flde Oliveira-Costa, 
Tegmark, Zaldarriaga, fc Ham ilton 2004 ; Sch warz, Stark- 
man, Huterer, fc Copi||20041 |Larid fc Magueijo l |2005| 
Copi, Huterer, Schwarz, fc Starkman | |2006[) , and trie low 
power at large scales (SpcrgeTet al.|2003 Copi, Huterer, 



Schwarz, fc Starkman|2U07| ); each oi these anomalies are 
<C 1% likely. For reviews with differing p oints of view 

see ~ 



Bennett et al. 



A natural quest 



|2010j |; |Copi et al.| ( pOlOl ). 
ion that arises is whether the align- 
ments and low power at large scales are correlated. 
Naively, the answer is negative, since the alignments are 
defined by orientation of the multipoles, and are indepen- 
dent of the total power. In the language of multipole vec- 
tors, the normalization at each multipole defines power, 
and is independent of the multipole vectors that define 
orientations of multipoles. Such a conclusion has been 



verified, for full sky maps, by Rakic & Schwarz (2007). 

However, the answer for the cut-sky maps is less clear 
a priori. Consider the case where the large-scale mul- 
tipoles are largely planar (as is the case for our own 
quadrupole and octopole), thus creating alignments. If 
the sky cut happens to be parallel to this plane, then the 



cut-sky power may be unusually low, and this may well 
be happening more often than when the cut is applied 
to unaligned (statistically isotropic) maps. Clearly, the 
only safe way to investigate the correlation between the 
low power and alignments in the presence of the realis- 
tic sky cut is to perform Monte Carlo comparisons with 
Gaussian random, statistically isotropic maps. 

In this Letter we investigate the correlation between 
the alignments and low power for the cut-sky as well as 
full-sky maps assuming a Gaussian random, statistically 
isotropic, cosmological model. 

2. LARGE-SCALE ANOMALIES AND WMAP7 

2.1. Low Power on Large Angular Scales 

The temperature anisotropy in the e [= (9, if))} direc- 
tion of the microwave sky can be represented by a real 
scalar function {AT(9, 4>)) on a sphere, expanded in 
terms of the multipole moments 



AT(e) = AT (6, 



£=0 



?HM), 



(i) 



where the £th multipole, 7| (#,</>), is expressed in terms 
of the spherical harmonics, Yg m (0, 4>), as 



Tr 



E 



(2) 



The standard inflationary cosmological model predicts 
that the fluctuations in the microwave sky can be thought 
of as being sampled from a statistically isotropic, Gaus- 
sian random field of zero mean. Gaussianity dictates 
that the variances of these ag rn would fully character- 
ize the distribution and statistical isotropy implies that 
these variances would depend only on £, allowing us 
to write the ensemble average of any pair of ai m as 



2 



(a£ m a* e , m ,) = Ci5u>5 mm > , where Ce is the (ensemble av- 
erage of) power in the £th multipole. In practice, it is 
not possible to measure such an ensemble average as we 
have only one realization of the universe. Instead, one 
uses its observable estimator, Ce, defined as 
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(3) 



In the case of statistical isotropy, this Ce, known as 
the angular power spectrum, is an unbiased estimator of 
Cg. Additionally, if one assumes Gaussianity, Ce is the 
best estimator of Ce with cosmic variance: Var(C^) = 
26$/ (21+ 1). 

Instead of calculating the power in each multipole, 
one can consider the (real space) two-point angular 
correlation function of the CMB temperature fluctua- 
tions. Assuming statistical isotropy, one can express 
the angular correlation function as C(ei,e 2 ) = C(9) = 
(T(ei)T(e 2 ))e, where (-)g represents an ensemble aver- 
age over the temperatures T(ei) and T(e 2 ) in all pairs 
of directions ei and e 2 separated by the angle 0. One 
can get an unbiased estimator of C(9) by replacing the 
ensemble average with the sky average 



C{9) = T(e x )T(e 2 ) 6 



(4) 



Note that Ce and C(9) only contain precisely the same 
information for full-sky data, and their analogues in the 
ensemble are informationally equivalent only if the sky is 
statistically isotropic. Thus, simultaneous measurements 
of both Cg and C(9) can be used to probe the validity of 
the assumption of statistical isotropy. 

The two-point angular correlation function, C(9), can 
be directly measured as there is a very large number 
of independently measured pixels on the WMAP sky 
In Fig. [l] we show C(9) = T(e-| )T(e 2 ) p (computed in 
pixel space, using SpICE (Szapudi, Prunet, & Colombi 
21X31) at NSIDE=512) for the WMAP 7-year coadded 
data. The results are shown for four different maps of 
WMAP7: the V and W bands masked with the 7-year 
KQ75 mask (KQ75y7; henceforth), the ILC map (which 
covers the full sky), and the cut-sky version of the ILC 
map (using the same mask). For comparison, we also 
show the full-sky ILC from the 5-year data. Finally, we 
show the angular two-point correlation function for the 
best-fit ACDM model for WMAP 7 data, along with the 
ler cosmic variance band (in blue) around the best-fit. 
The most striking feature of C(0) for all the cut-sky maps 
is that they are very close to zero for 9 > 60° , except for 
some anti-correlation near 180°. 

To quantify this lack of correlation, we follow past work 
dSpergel et all |2003| [ Cggi et al.||2007[ |Copi, Huterer7 
Schwarz, & Starkman||2l)U9|) and adopt the Si/ 2 statistic 
to test the total amount of correlation at angles above 
60°: 



Sl/2 = 



1/2 



C(0) d(cos0). 



(5) 



The calculation of S1/2 using Eq. (j5j) is susceptible to 
small-scale fluctuations in C(0). To avoid this we evaluate 
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Fig. 1. — The two point angular correlation function from the 
WMAP 7-year coadded data. The long-dashed and the dot-dashed 
lines show C(9) = T(ei)T(e2)g for the V and W bands, respec- 
tively, both masked with the KQ75y7 mask. The dashed line shows 
the correlation function for the ILC7 map (which covers the full 
sky). The solid line shows the same for ILC5 for comparison. The 
C(0) for the cut-sky version of the ILC7 map (using the same mask) 
is shown by the dot-dot-dot-dashed line. Finally, we show the an- 
gular two-point correlation function for the best-fit ACDM model 
with the dotted line, the ltr cosmic variance band around the best- 
fit being shown in blue. 



S 1/2 via ( |Copi et al.|2009[ ) 

Sl/2 



JL^Jtf+lXJtf + \)CeIe.e- (1/2)4, (6) 



where h,i'{x) — j_ x Pe.(x')Pii(x')dx' and Ce is the Leg- 
endre transform of C{9) 



Ce = 2tt 



P^(cos0)C(0)d(cos0). 



(7) 



The S1/2 values, calculated using Eq. (j6|, for different 
bands of WMAP 7-year (as well as 5-year, for compar- 
ison) coadded data is shown in Table [ll Here we make 
use of SpICE at NSIDE=64. 

2.2. The Quadrupole-Octopole Alignment 

To measure the level of alignments between the 
quadrupole and octopole, we use two of the most com- 
monly used statistics: one based on the multipole vec- 
tors, and one based on the angular momentum opera- 
tor. As expected, the two statistics are mutually highly 
correlated (see Table [2] below) ; we include both for com- 
pleteness. 

Instead of expanding the temperature anisotropy of the 
microwave sky in spherical harmonics, one can uniquely 



direction as (Copi, Huterer, & Starkman||2004 


T t (0,<j>)=AW 







(8) 



where is a scalar which depends only on the total 
power in the I th multipole and {v^'^|z = 1, ...,£} are I 
'headless' unit vectors, known as the "multipole vectors" . 
Here Te is the sum of all possible traces of the first term 
in the right hand side of the Eq. 
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TABLE 1 

-STATISTIC FROM WMAP 5-YEAR AND WMAP 7- YEAR DATA 



Data 
Source 


V5 
(KQ75y5) 


W5 
(KQ75y5) 


ILC5 
(KQ75y5) 


IT C5 V7 
ILC5 (KQ75y7) 


W7 
(KQ75y7) 


ILC7 

(KQ75y7) 


ILC7 Theory5 C e 


Theory7 C e 


Sl/2 

(MK 4 ) 


1242 


1236 


1054 


8600 1312 


1231 


1141 


8484 49040 


46610 



Note: The values of theSj/ 2 -statistic are calculated for different bands of WMAP 5-year and WMAP 7-year coadded data. We have used 
SpICE | |Szapudi et al.||2001| at NSIDE=64 and calculated S 1 / 2 using Eq. j6|l with £ m ax = 30. The values of S1/2 calculated from the 
best-fit theory Ui are also included for reference. 



To test the planarity of the quadrupole and octopole, 
for each multipole £ we form t he £(£— 1)/2 cross-pr oducts 
(the "oriented area" vectors ( Schwarz et al.| 2004)) 

w Wj)=*m x *<M) t (9) 

where the overall signs of the area vectors, w^ ;lJ \ are 
unimportant. For the following discussion, let us con- 
sider the lone area vector for the quadrupole, w( 2;1 ' 2 ), 
and the three area vectors for the octopole: w^ 3;1 ' 2 ', 
w^-^andw^ 1 ). 

To investigate the relative orientation of the quadrupo- 
lar plane with the three octopolar planes, we can evaluate 
the magnitudes of the dot products between w^ 2 ' 1,2 ) and 
w (3;»j) j given by 



A\ {or 2 or 3} 



A 2 ' 1 ' 2 ) . w( 3 ^' 2 ) f° r ( 3 ; 2 ,3) or ( 3 ;3,i)} 



(10) 

High values of A4 imply that the £ = 2 and £ = 3 planes 
are aligned with each other. One can test the combined 
planarity of the quadrupole and octopole by using the 
statistic that takes the average of the dot products 

1 



5(3,9) 
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Using ILC7 and ILC5, we find = 0.736 and S( 3 > 3 ) = 

0.753, respectively. In the subseque nt sections , to be 



conservative, we use 5( 3 > 3 ' = 0.798 dCopi et al 



Legm 



, 20071), 

ark. 



obtained using the cleaned full sky map ot [T 
|Qliveira-Costa, & Hamilton| ( |2003f (TOH; henceforth) 



c e 



An alternative statistic that tests th e planarity is given 
by th e angu lar momentum dispersion (de Oliveira-Costa 
et al. (2004)). For each £, an axis can be found, around 



which the angular momentum dispersion 

1 

(T e (h e )\(h e -Lf\T e (h i ))= J2 m 2 \a £m (n e )\ 2 (12) 

m=-t 

is maximized. Here L is the angular momentum operator 
and ai m {hi) are the spherical harmonic coefficients of the 
CMB map in coordinate system with its z-axis in the fig 
direction. Note that this angular momentum axis distills 
a limited amount of information from each multipole, 
r educing the d.o.f. t o just 2 from the 2£ d.o.f. in v^'^ 



Copi et al (|2006j) found that n 2 • n 3 = 0.962 for the 
TOH map. The corresponding values obtained from the 
analysis of the WMAP 7-year data is 0.937. The mea- 
sured values of S( 3 ' 3 ) and A2 • A3 consistently indicate 
that the quadrupole and the octopole are unexpectedly 
aligned, both in the 5 and 7 year WMAP data. 

3. IN SEARCH FOR A CORRELATION 



We have just seen that WMAP7 confirms results from 
its earlier data releases: at large angular scales (> 60°), 
the two-point correlation function of the temperature of 
the CMB, as measured by WMAP, is significantly smaller 
in magnitude than expected from the predictions of the 
standard inflationary cosmological model; and in addi- 
tion, the planes of the quadrupole (I = 2) and the oc- 
topole (£ — 3) are unexpectedly aligned with each other. 
In this section, we investigate whether these anomalies 
are correlated at a statistically significant level. 

More specifically, we want to investigate a couple of 
different scenarios: 

• If the Gaussian random, statistically isotropic 
CMB maps are constrained to have low angular 
correlation on large angular scales, are they more 
likely to exhibit planarity and alignment of the 
quadrupole and octopole? 

• If the Gaussian random, statistically isotropic 
CMB maps are constrained to have aligned 
quadrupole and octopole, are they more likely 
to have low angular correlation on large angular 
scales? 

Although the answers to the above questions might 
seem to be obvious for the full-sky case, there is no direct 
reason why no such correlation should exist for the case 
of cut-sky maps. In what follows, we attempt to address 
each of these questions for both the full-sky and the cut- 
sky cases. 

3.1. Monte-Carlo Analysis 

We first generate 100,000 constrained (low angular cor- 
relation on large angular scales) realizations of the CMB 
sky with selection criterion of S1/2 < 8583(/x.?f) 4 (i.e 



lower than the ILC5 full-sky S 1/2 from Copi et al. ( 2009 )) 
and the cut-sky S1/2 < 1 1 52 (/xi'T) 4 (i.e. lower than the 
ILC5 cut-sky map)F] This is our sample of Gaussian ran- 
dom, statistically isotropic maps that are constrained to 
have low angular correlation on large scales, and we use 
them to test the probability of alignments under this con- 
straint, i.e., the probability of having 5^ 3,3 -' > 0.798 for 
these low-correlation maps. We also calculate the proba- 
bility of these low-correlation maps to have aligned angu- 
lar momentum dispersion axes, i.e., n 2 • 113 > 0.962. We 

It should be noted that these values of S1/2 are corresponding 
to the WMAP 5-year ILC full sky (DQ corrected) and WMAP 5- 
year ILC cut sky (using KQ75y5 mask, DQ corrected) maps, since 
we started this analysis prior to WMAP 7-year data release. In 
fact, these values are comparable to the corresponding values for 
the WMAP 7- year results shown in Table [l] 



4 



TABLE 2 

Probabilities of having Low-Power on Large Angular 
Scales and Quadrupole-Ootopole Alignment for the 
Unconstrained and Constrained Maps 



Probabilities 
(in %) 



PlS 



full sky 
1/2 



-rcut sky 
'l/2 
p'( S (3,3)) 

P (n2 ■ A3) 




Constrained Maps 
Low Power Full-Sky Aligned 



0.12 
0.36 



6.9 
0.07 

99.6 



Note: The statistics M 3,3 ) and n2 -113 are always evaluated on the 
full sky. The cases with the dashes denote circular comparisons. 

then compare these probabilities with the corresponding 
ones for the unconstrained maps. 

To test the probability of low angular correlation on 
Gaussian random, statistically isotropic maps that are 
constrained to be aligned, we generate 100,000 MC maps 
constrained to have S 1 ' 3,3 ' > 0.798, calculated using the 
full sky as the cut sky alignments are extremely dif- 
ficult to test (the errors in multipole vectors become 
large for a cut larger than a few degrees). We then 
evaluate the Si/ 2 statistic for these constrained maps 
(with and without applying the WMAP 5-year KQ75 
mask (KQ75y5)) and calculate the percentage of these 
maps having Si/ 2 < 8583(/i-ftT) 4 (for the full-sky case) 
or S 1/2 < 1152(Vo 4 (for the case using KQ75y5). We 
also calculate the probability of these constrained maps 
to have n 2 • n 3 > 0.962. As discussed before, we then 
compare these probabilities with the corresponding ones 
for the unconstrained maps. 

3.2. Results and Discussions 

Our main results are summarized in Tabic [2] Using our 
sample of 100,000 low-power maps, we find the probabil- 
ity of these maps to have S^ 3 ' 3) > 0.798 is 0.12%, which is 
the same as the probability of completely un constrained 
maps to have quadrupole-octopole alignment ( Copi et al. 
2006 1 . This shows that the low-power maps do not have 
a higher probability to have an alignment of the I — 2 
and i = 3 planes. The probability of these low-power 
maps to have aligned angular momentum dispersion axes 
(n 2 • n3 > 0.962) is also very low (0.36%) and essentially 
identical as that found for unconstrained maps (0.37%; 



Copi et al. (|2006 )). 

Using 100,000 maps constrained to be have 
quadrupole-octopole alignment (S^ 3,3 - 1 > 0.798), 
we find 6.9% probability for these maps to have 
"5*1/2 < 8583(/iif) 4 for the full-sky case. Comparing 



this with the probability (7.0%) of having low angular 
correlation for completely unconstrained maps, we 
find no correlation between alignment and low angular 
correlation for the full-sky case. 

For the cut-sky case (using the KQ75 mask), how- 
ever, we find 0.07% probability for these maps to have 
Si/2 < 1 152 (/iif) 4 . This probability is slightly higher 

(> 2cr with a a Poisson error) than P {s^l^ = 0.05% 

for unconstrained maps. However, this difference is not 
statistically significant to conclude any detectable corre- 
lation between alignment and low power for cut skies. 
4. CONCLUSION 

Applying the multipole vector formalism to Gaussian 
random and statistically isotropic realizations of CMB 
maps (both full-sky and cut-sky) lacking correlation on 
large angular scales, we find no increased probability for 
alignment of the quadrupole and the octopole than ex- 
pected in the case of unconstrained random maps. These 
low-power maps also do not show higher angular mo- 
mentum dispersion. On the flip side, we also find that 
realizations of CMB maps (both full-sky and cut-sky), 
constrained to have aligned quadrupole and octopole at 
a level greater than that exhibited by WMAP, do not 
have lower power on large scales than expected in the 
case of unconstrained random realizations. 

Our results conclusively demonstrate that, under the 
standard Gaussian and isotropic model, there is no statis- 
tically significant correlation between the missing power 
on large angular scales in the CMB and the alignment 
of the 1 — 2 and £ — 3 multipoles. Therefore, in the 
context of the standard model, their combined statisti- 
cal significance is equal to the product of their individ- 
ual significances. For example, simultaneous observation 
of the missing large-angle correlations with probability 
P(S'™2 Sky ) ^ 0.1% and alignments with the probabil- 
ity P(S*( 3 ' 3 )) ~ 0.1% is likely at the < 0.0001% level. 
Given that both anomalies occur at the largest observ- 
able scales and are correlated with special directions in 
the sky (ecliptic and/or dipole), they clearly require a 
causal explanation. 
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